UG-DFM-2/97 
FTUV/97-18 
IFIC/97-18 



Many Body approach to the inclusive (e, e') reaction 
from the quasielastic to the A excitation region 



A. Gil\ J. Nieves^ and E. Oset^ 

^ Departamento de Fisica Teorica and IFIC, Centra Mixta Universidad de Valencia 
CSIC, 46100 Burjassat (Valencia) Spain. 

^ Departamento de Fisica Maderna, Universidad de Granada, 18071 Granada, Spain. 



Abstract 

We have performed a many body calculation of the inclusive (e, e') cross sec- 
tion which runs over the three traditional regions at intermediate energies: the 
quasielastic peak, the dip region and the delta region. The longitudinal and 
transverse response functions in the quasielastic peak have also been evaluated. 
Traditional effects like polarization, meson exchange currents, final state interac- 
tion and delta renormalization in the nuclear medium have been included. Meson 
exchange currents are generated from a model of pion electroproduction on the 
nucleon which reproduces accurately the experimental data. 

The inclusive cross section accounts for lA^, 2N, 3N mechanisms of virtual 
photon absorption and one pion production. Meson exchange currents associated 
to the (7*,27r) reaction are also accounted for. 

We obtain good results for the (e, e') cross sections in the whole energy range 
and for different nuclei. The response functions are also in good agreement with 
the latest experimental analysis. On the other hand, the method provides the 
separation of the contribution to the inclusive cross section from different phys- 
ical channels which is a necessary input to evaluate cross sections like {e,e'N), 
{e,e'NN), (e.eV) etc. 



1 Introduction 

Inclusive electron scattering, particularly around the quasielastic peak is probably one 
of the problems that has attracted more attention in nuclear physics. One of the 
reasons for the proliferation of theoretical work was the persistent difficulty of simple 
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pictures based on the shell model of the nucleus, or the Fermi gas approach, to repro- 
duce simultaneously the longitudinal and the transverse response functions around the 
quasielastic peak. Particularly, the longitudinal response appeared systematically to 
be much larger than experiment [|I], ^ 

Even more worrisome was the fact that the integrated strength of the experimental 
longitudinal response function showed large discrepancies with the expected result 
according to the Coulomb sum rule [Q. At large values of the momentum transfer q 
this integral should be the charge of the nucleus and there was some apparent missing 
strength. 

Those experimental results have been generally accepted, even when some other 
experimental results seemed to challenge them. Indeed the results of the experiment of 
1^ in 238u j-^q|- ^Y^fy^ ^Y^Q expected suppression of the longitudinal response. More 
recently an experiment at Bates on ^''Ca [^] showed a longitudinal response larger than 
in ref. [Q with only 20% reduction over the simple shell model expectations. 

The growing experimental discrepancies stimulated the thorough and thoughtful 
work of [0. In this work the author analysed the world set of data and made further 
improvements in the analysis, coming with new results for the response functions which 
show a smaller reduction of the longitudinal response than previously assumed. At the 
heart of the issue was the fact that, in a modified Rosenbluth plot in terms of the 
variable e (which runs from to 1), the Saclay data ||, Q concentrated their points 
in the region of e < 0.5 which induces large errors in the slope of the straight line which 
correlates the points of the plot. Those results, complemented by others at SLAC and 
Bates, which fill up the region of e ~ 1 , lead to a more accurate determination of the 
slope and thus the longitudinal response. The problems with the Coulomb sum rule 
are then automatically solved to the relief of all p. 

As usually happens, some theoretical calculations were "successful" in explaining 
the abnormal reductions of the longitudinal response, through the use of small effective 
masses, swelling of nucleons, abnormal nucleon form factors, correlations, etc. But in 
the benefit of the theoretical community it should be said that whenever this happened, 
the agreement with the transverse response was spoiled, although excuses and good 
wishes beyond the ideas used for the longitudinal response were invoked as possible 
solutions. Hence, a convincing simultaneous explanation of both the longitudinal and 
transverse response functions was never obtained. 

The interesting thing is that the existence of a "problem" induced so much work 
that practically all resources of nuclear physics and many body theory have been used 
in this topic and a lot of things have been learned. In the present work we shall benefit 
from all this previous work and by using a selfcontained many body formalism we 
will incorporate all these effects which convincingly proved to be relevant in previous 
works. We also incorporate some new ingredients which come out naturally within 
our many body expansion and furthermore we include Delta-hole (A/i) excitations, 
additional to the particle-hole {ph) excitations, which allows us to simultaneously 
study the quasielastic peak, the A peak and the "dip" region between the two peaks 
in the inclusive (e, e') cross section. 
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2 Brief review of different approaches and intro- 
duction to our approach. 

The amount of ideas which have been studied in connection with the inclusive (e, e') 
scattering is quite large. We shall discuss them briefly: 

2.1 Modification of the nucleon form factor (swelling of the 
nucleon) 

These ideas were soon invoked and explored within different frameworks 0, 1^. They 
appear naturally in some microscopic models of the nucleon form factor when the 
intermediate baryon propagators are replaced by those in the nuclear medium. This is 



the case, for instance in ref. [|TT|, where the underlying elementary model is the Nambu- 
Jona-Lasinio model of the nucleon. Obviously there are many other medium effects 
not taken into account in these schemes, as we shall see. Furthermore, although in a 
different language and in terms of a few relevant physical magnitudes, such medium 
modifications appear in a systematic many body expansion, where they can be classified 
as vertex corrections. 



2.2 Relativistic effects 



Theories using relativistic scalar and vector potentials like in the Walecka model ||T2 



have been popular. The scalar and vector potentials are about one order of magnitude 
larger than the ordinary non relativistic potential which is roughly the sum of the two. 
This cancellation is missed in many applications of these relativistic potentials leading 
to unrealistic predictions. The appealing thing of the relativistic approach is the small 
effective mass, ~ Mj^ j2 of the nucleon and the fact that the nucleon response for 
■ph excitations is roughly proportional to M]^. This reduces the longitudinal response 
but also the transverse one. 

A clarifying view of these problems is exposed in where the necessity to go 
beyond the relativistic mean field approach is shown in order to avoid the pathological 
predictions tied to the small effective masses, like in the computation of the nuclear 
response functions or the large relativistic enhancements that drive magnetic moments 
outside their Schmidt lines. In |TB] a relativistic RPA calculation is used. Similar con- 
clusions are found in stating that "selfconsistent calculations" show cancellations 
between large relativistic effects on the single nucleon current and on the many nucleon 
wave functions. In [jl5[ it is also shown, by solving numerically the Dirac equation with 
the relativistic potential, that the genuine relativistic effects in the enhancement of the 
axial charge amount to 20-30%, while perturbative calculations give as much as 70-80% 
enhancement 
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Even when improved with some selfconsistent steps, relativistic calculations still 

which is only an approximation to 



rely on the concept of the effective mass |1T3| , p!8| , |T9 



the richer content of the nucleon self-energy. The nucleon self-energy is a function of 
the energy and momentum, as independent variables, and leads to important dynam- 
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ical properties of the nucleus [20| not contained in static pictures like the mean field 



theories. One of the consequences is that the nucleon effective mass is a strongly depen- 



dent function of the energy, with a peak around the Fermi surface [20|. Furthermore 



as shown in |^T|, RPA correlations tied to the pionic degrees of freedom are essential 
in order to provide the energy dependence of the nucleon self-energy and hence the 
dynamical properties of the nucleus. Actually it is quite interesting to see that recent 
sophisticated calculations in light nuclei using path integral Monte Carlo methods p2| 



find that "pion degrees of freedom in both nuclear interaction and currents play a 
crucial role in reproducing the experimental data" . 

From this discussion it looks clear that improvements along the relativistic line for 
the present problem should include "selfconsistency" in the sense of ref. |T^, in order 
to exploit the large cancellations between the scalar and vector potentials. 

2.3 Pionic effects 



As mentioned above the pionic degrees of freedom play an important role in 
quasielastic electron scattering. This has also been emphasized in where meson 



exchange currents driven by pion exchange are evaluated, putting special emphasis in 
fulfilling the continuity equation and preserving gauge invariance in the many body 
system. This imposition has as a consequence some changes in the results with respect 



to former works along similar lines [|^, |2^, |26| 



In ref. ||2^ similar ideas, but using the formalism of path integrals is followed. RPA 
correlations are automatically generated in that scheme leading to some quenching of 
the longitudinal response from a reduction in the isoscalar channel. 

Pions are also explicitly used in approaches which include meson exchange currents, 
as we shall see below. They are usually taken static, as in meaning that the 

energy carried by the pion is neglected. While this is a fair approximation for the 
exchange currents at energies below pion production threshold, at higher energies the 
need to work will the full pion propagator becomes apparent. This is particularly true 
if one wishes to account for real pion production in the same many body scheme, as 
we shall do. 

In the resonance region primary pion production accounts for the largest part of 
the response function (although some of the pions are absorbed in their way out of the 
nucleus and show up in 2N or 3A^ emission channels). Hence, the explicit treatment 
of pionic degrees of freedom allowing pions to be produced, both as virtual as well as 
real states, becomes a necessity in this region. Our scheme puts a special emphasis 
on pions. In fact it follows a different path to other schemes, beginning with real 
pion production and ensuring that a proper hand on the (e, eV) reaction is held. Then 
exchange currents and further corrections in the many body system are generated from 
the model for the eN — > e'A^yr reaction. 

2.4 Meson exchange currents 

A large fraction of work has been devoted to the role played by meson exchange cur- 



rents (MEG) in this reaction |2^, ^ |2^, ^ 0. The standard seagull, pion 



4 



in flight and A terms are included mediated by pions. The indirect effect of short 
range correlations in these terms is, however, neglected in those works. The work of 
3T| incorporates terms in the scheme which account for virtual photon absorption on 



correlated pairs, hence accounting for ground state correlations. Although the same 
concepts are shared in the previous approaches, differences in the input and the way 
to implement them lead to different results. 

Our approach differs from the quoted works although conceptually it is quite similar. 
First we realize that the two body currents appear in (e, e') as corrections to the main 
one body contribution. This is because virtual photons can be absorbed by one nucleon. 
This is opposite to the case of real photons which require at least two nucleons to be 
absorbed (we are thinking in terms of infinite nuclear matter). Thus it is clear that the 
laboratory to test the effect of two nucleon currents is real photons not virtual ones. 
This is the reason why prior to the present work we devoted energies to the problem of 
real photon absorption |^ . Second, in order to minimize sources of uncertainties, the 



MEG were generated from the model nN by allowing the pion to be produced 

in a virtual state and be absorbed by a second nucleon. After this is done, long range 
correlations from polarization phenomena, as well as short range correlations, are taken 
into account. The model for the — >■ t^N reactions was tested against experimental 
data and was found to be good. This gives one some confidence in the strength that 
one generates for the MEG. The reliability of the method gets extra support from 
recent measurements of two body photon absorption |3^, |3^, ^ where the agreement 
with the (7, np) emission channels is rather good. Some discrepancies remain in the 
{liPP) channel, but this channel has an experimental cross section nearly one order 
of magnitude smaller than the (7, np) one, thus for the purpose of the total (7, NN) 
emission or the two body MEG in (e, e') such discrepancies will not play an important 
role. 

In view of the success in the real photon case, we adopt here the same scheme and 
follow the same steps simply substituting the real photon by the virtual one. This 
means that we begin with a model for the eN — >■ e'Nir reaction and construct the 
MEG from it following identical steps to those of ref. For this reason we begin 



in next section by showing our model for the eA^ e'Nn reaction and contrasting it 
with the experimental results. 

2.5 RPA correlations 

Several works have emphasized the role played by RPA correlations allowing for a 
ph excitation which propagates in the nuclear medium mediated by some residual ph 



interaction |3^, ^ . These works share the feature that a reduction is produced 

in the longitudinal channel. We shall also incorporate these long range correlations or 
polarization effects. In addition we shall also include Ah excitations, as a source for 



polarization which will be relevant in the transverse channel as shown in [38 
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2.6 Final state interaction (FSI) 



This is another topic which has received some attention and a thorough work devoted 
to this issue can be seen in It is clear that once a ph excitation is produced by the 
virtual photon, the outgoing nucleon can collide many times, thus inducing the emission 
of other nucleons. A distorted wave approximation with an optical (complex) nucleon 
nucleus potential would remove all these events. However, if we want to evaluate the 
inclusive (e, e') cross section these events should be kept and one must sum over all 
open final state channels This is done explicitly in and the result of it is 

a certain quenching of the quasielastic peak of the simple ph excitation calculation and 
a spreading of the strength to the sides of it, or widening of the peak. The integrated 
strength over energies is not much affected though. 

The use of correlated wave functions, evaluated from realistic A^A^ forces and incor- 
porating the effects of the nucleon force in the nucleon pairs has also been advocated 
in connection with the effects from final state interaction ETl . If one incorporates two 



particle-two hole (2p2h) components in the final excited states one gets the spreading 
of the peaks as found in . For the purpose of the response function it is like exciting 
ph components which have a decay width into the 2p2h channel. This gives rise to the 
quenching of the peak and spreading of the strength. Another source of these effects 
is the momentum dependence of the nucleon self-energy which is also accounted for 
in the scheme and which sometimes is taken into account approximately in terms of 
an effective mass (although our position on evaluations which use this variable has 
been already made clear). The approach of [JT| using an orthogonal correlated basis 
with functions obtained with variational methods also incorporates RPA correlations 
discussed in point 5) but does not account for the MEC discussed in point 4), which 
rely upon the coupling of the photon to the pion or the A excitations. 

In our many body scheme we will account for this FSI by using nucleon propagators 
properly dressed with a realistic self-energy in the medium, which depends explicitly on 
the energy and the momentum This self-energy leads to nucleon spectral functions 
in good agreement with other accurate more microscopic approaches like the ones in 
The self-energy of has the proper energy-momentum dependence plus an 



44 



imaginary part from the coupling to the 2p2h components, hence it has the ingredients 
to account for the FSI effects discussed in although using a different calculational 
scheme. 

Nuclear spectral functions with a language closer to the one we shall follow have 
also been used in They include the interesting result, which we will employ here, 
that keeping the width of the particle states is important but one can disregard the 
width of the hole states. 



2.7 A excitation 

While many efforts have been devoted to the quasielastic peak very little attention has 
been given to the A region and the dip region between the quasielastic and A peaks. 
One exception is the work of |^6[ which looks at the effects of MEC in the dip region. 
A recent work presents some experimental results and a theoretical analysis of 
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the A region based on the eA^ —>■ e'Nir model of |4^. One interesting conclusion of the 
work is that the data have a broader energy spectrum than the theoretical one based 
on the properties of a free A width. This suggest a larger A width from A coupling 
to many body channels, additional to the natural decay width with effects of Fermi 
motion included. This is actually well known from pion physics |]49| , |50| . 

Our aim in this section is also to evaluate as accurately as possible the response 
function in this region, for which we shall use results for the A in a nuclear medium 
^T|] which have been tested thoroughly in a variety of pionic reactions: elastic 



quasielastic, charge exchange, absorption, etc... [§3 . 

The discussion in this sections has served to expose relevant works done in the 
literature and to present our approach in connection to the ideas exposed in these 
works. We shall follow a microscopic many body description of the inclusive (e, e') 
process and will incorporate in our approach the important effects discussed in this 
section. We are thus aiming at an evaluation as accurate as possible at the present 
moment. A small sacrifice is made in order to allow one to treat microscopically in 
a tractable way all the effects discussed above: we use an infinite nuclear medium 
and obtain results for finite nuclei using the local density approximation (LDA). This 
was also used for the evaluation of the total cross section of real photons with nuclei 
with good results |^ and a mathematical derivation was made there to justify the 
accuracy of the approximation. A direct comparison of results with the LDA with 
finite nucleus results, using the same input was done in for deep inelastic lepton 
scattering, showing that indeed the LDA is an excellent approximation to deal with 
volume processes (i.e no screening or absorption effects). In a comparison of the 
results for finite nuclei with those of the Fermi gas around the quasielastic (e, e') peak 
is done, showing that for some average Fermi momentum, the results of the Fermi gas 
and those of finite nuclei are nearly identical, proposing that choice as an even better 
prescription than the LDA. We shall follow the LDA, already tested for real photons, 
and sufficiently good for our purposes. 

Another feature which is novel in our approach and rather important in practical 
terms is the following: we shall use a method which relates the (e, e') cross section 
to the imaginary part of the virtual photon self-energy. This is actually not new 
and has been used before [0, The novelty is that by using properly Cutkosky 
rules one can relate the different sources of imaginary part to different channels which 
contribute to the inclusive (e, e') cross section. In this way, we lay the grounds to 
evaluate from the present input the cross sections for exclusive processes {e,e'N), 
{e,e'NN), (e, e'vr) (e, e'A^vr) etc., which will be the subject of a forthcoming paper 



The results of the present paper and those of are part of a PhD thesis [pS 
where additional details to those given here can be found if desired. 



3 The eN e'Nn reaction 
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3.1 Formalism 



We shall follow the model used in ref. |^ for ttN at intermediate energies 

generalizing it to virtual photons. In fact this model is essentially the same as used 
for the eA^ — eNn reaction in . However, having in mind the application to nuclei 



we make a reduction of the relativistic amplitudes keeping terms up to with, 
P,Miy the nucleon momentum and mass. The neglect of the terms is justified 

numerically as we shall see in the results, so we construct this non relativistic amplitude 
for the eA^ — > eNn process ready to be used with ordinary non relativistic nucleon wave 
functions. 




N 



N 



(b) 




Fig. 3.1 Basic couplings of the virtual photon and the pion to the nucleon and to the A^A 
transition. 

The basic couplings which we need are those depicted in fig. 3.1 which account for 
the coupling of the photon and the pion to the nucleon and to the A^A transition, plus 
the KroU Ruderman term (KR), and the coupling of the photon to the pion. The KroU 
Ruderman term appears as a gauge invariant term through minimal substitution when 
a pseudovector ttA^A^ coupling is used, as we do. The analytical expressions for these 
vertices are given in appendix A. For convenience, from the KR term of the appendix 
coming from gauge invariance, and which we call there seagull term, we construct the 
KR term displayed in this chapter such that it contains all the non vanishing pieces of 
the amplitude when — > 0. We will come back to this point later on. 



3.2 The amplitudes for the eA^ eNn model 



The Feynman diagrams which are considered in the model for 7A^ irN of [32| or in 
7*A^ — s> vrA^ of [^] (7* will stand from now on for the virtual photon) are depicted in 
fig. 3.2. 



8 



N P 




(a) 




(b) 
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(d) 



(e) 



(0 



Fig. 3. 2 Feynman diagrams considered for the 7*A^ ttN reaction. 

They are the nucleoli pole direct (NP) term (a), the nucleon pole crossed (NPC) 
term (b), the pion pole (PP) term (c), the delta pole direct (DP) term (d), delta pole 
crossed (DPC) term (e) and Kroll Ruderman (KR) term (f). The expressions for these 
amplitudes are obtained by doing the nonrelativistic reduction of the relativistic ampli- 
tudes of . There is some small contribution from the negative energy intermediate 
nucleon states which is kept in our expressions, which as we mentioned neglect only 
terms of order O(-tj-)^. The corresponding expressions are: 
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X 



FN'( 2^ j 2p + q-2k \ o {axq)ak 



with 



If we consider, 



we get 



B{n,n7r^) = -1 
B{n,pTi-) = V2 

B{p,n7r+) = V2 



/(TfO) = /.(TfO) = 2/3 
/(7r+) = -/,(7r+) = -v^/3 
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where ^/s is the invariant energy of the system virtual photon-initial nucleon, e(e > 0) 
is the electron charge, m^r = 139.5 MeV is the pion mass, Ma = 1238 MeV, the A 
mass 



32| and the free decay width of the A is given by 



1 / f* \ ^ M 

r(^) = 7r — ^\kcm\'QiV^s-Mr,-m^) (7) 

On the other hand, /t^nn, f.y, f* , Fi,Gm, Fa, F^^.^,, Ft,, are coupling constants and 
form factors which we show in the appendix, as well as S, the spin transition operator 
from spin 3/2 to 1/2. 

The expressions given keep the Lorentz covariance up to terms The small 

A crossed term of eq.(6) holds strictly in the frame where the outgoing nucleon is at 
rest. 

As mentioned before, in the KR term of eq. ( ^ we have included corrections in the 
zero component which come from the zero component of the NP and NOP terms when 
using the vertex of eq.( |104|) (this means from positive energy intermediate states). 
On the other hand the q^/2MN terms in the spatial components of the KR term come 
from the intermediate negative energy components of the NP and NCP relativistic 
amplitudes. This trick serves us to concentrate on the KR term all the contributions 
which do not vanish when the pion momentum goes to zero. 



3.3 Gauge invar iance and form factors 

Gauge invariance implies that q^M.^ = where Ai'^ = J^i-^i- From the expressions 
given above we can see that the delta terms are gauge invariant by themselves. As 
for the rest of the amplitudes they form a block of gauge invariance terms in the 
absence of form factors. However, we must impose some restrictions in order to keep 
gauge invariance when the form factors are included. If we consider the 7*n — pir^ 
amplitude we see 



el^B{N,N'^) 



P 



r - E{p- k) 



(8) 



which together imply 



-aqFA - F^,,,,a{k - = 



Fl{q') = F,^^{q') = F^iq' 



(9) 



Gauge invariance in other isospin channels does not require extra relationships. In 
the Appendix we can see that A ~ Ma — V^Pn and hence eq. (9) is fulfilled to a good 
degree of approximation. However, in order to keep strict gauge invariance we take 
only one form factor for all, which we choose to be Ff of eq.( |108|) . We have checked 
that by taking any of the other form factors the changes induced in the cross sections 
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are much smaller than the experimental errors (see fig.3.6), so we are rather safe with 
any of these choices. 

On the other hand, since in the pion pole term we have a form factor corresponding 
to the coupling nNN with a virtual pion, — k)^), we have included this form 

factor in the amplitudes of this block (NP,NCP,PP,KR) in order to preserve gauge 
invariance. 



3.4 Unitarity 



Another refinement introduced in the model is unitarity. Watson's theorem implies 
that the phase of the ttN ttN and tcN — > amplitudes in each term of the 
partial wave decomposition must be the same. Such as our model stands, we have a 
A term which by itself satisfies the theorem if we assume the ttN — > vrA^ amplitude 
dominated by the A, as it is the case. However, in the ttN amplitude we 

have a sizeable background that in our model is real and the sum of the terms does 
not satisfy Watson's theorem. Although this violation of unitary does not result in 
important numerical changes in the cross section, we nevertheless unitarize the model 
as was done for real photons . We follow the procedure of ||5^ introducing a small 
phase (pi^/s, q^) which corrects the A term, where y/s is the invariant energy of the 
virtual photon-initial nucleon system and the four-momentum squared of the virtual 
photon. By means of an iterative method we find 0(v^, q^) such that 



Im 



(TAiq 



+ Ts{q')) 



(3/2,3/2) -i<5(3/2,3/2)(q2) 







(10) 



where T^iq'^) represents the A pole direct term amplitude, Ts^q^) is the contribution 
of the rest of the terms to the (3/2, 3/2) channel (see details in for the projection 
of these terms in the (3/2, 3/2) channel) and 5(3/2,3/2) (?^) are the vrA^ phase shifts in 
the 3/2,3/2 spin-isospin channel. 



3.5 Cross sections for the eN eNir process 



We follow here the steps and nomenclature of ref. In fig. 3. 3 we show diagram- 

matically the process with the different variables which we use, /Ce, fcg, q,p,p' and k 
representing the momenta of the incoming electron, outgoing electron, virtual photon, 
initial nucleon, final nucleon and pion respectively. 
The unpolarized cross section is given by 



r/^rr n"^ \k'\ 

^-VP\L,Ae,e')Wr^, (11) 



dn'JE', q' \l 
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where a = 1/137.036 is the fine structure constant (e^/47r) and L^j, the leptonic tensor 
defined as 

L^u{g, ^ ) ~ '^{kf.^kf.y + fcgj^fce^ + —Qfiu) (12) 
The hadronic tensor is given by 

^™ (27r)3 E' J {2nY2E^^''^ 

with j^^ the 7*A^ to A^'vr amphtude defined in sect. 3.2. 

We can separate from there the angular dependence of the pion and get 

"''^''-V(e,e')l^i!L(iV) (14) 



where now 



(27r)3 E' J {2'Kf2Ey^ 



(13) 



(15) 



Now by taking q along the z direction, using gauge invariance and the explicit 
expressions for L^^ we can write, following exactly the same steps as in ^ 
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(16) 

+72e(l + 6)^co.$;| 

where 



and 



r = 



a 




1 ■ 




27r^ 


ke 1 




l-e 


s — 









dar _ _e^_MiVit i \( jxx , jyy\-\CM 



(17) 



(18) 
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A;, 



■CM 



where the variables of the electron are in the lab frame while those of the pion are in 
the 7*A'" CM frame. J^'' in the former expressions is given by 



Ji^ — Tri^jJ^j^^ (19) 

The angular variables dg.-, 0*, $* are depicted in fig. 3.4. The variables ke, k'^ deter- 
mine the (e, e') reaction plane and q, the virtual photon momentum, determines the 
z direction. The q direction and the pion momentum k determine the irN reaction 
plane, and the angle between this plane and the (e, e') plane is the angle $* . 

The normalization of dai/dD,* is chosen in such a way that in the limit of real 
photons dai/dVl*^ coincides with the unpolarized cross section of 7iV — ^ -kN with real 
photons. In this case the variable k^ becomes the lab momentum of the real photon. 
The variables ctt, cr^, Up, o"/ are the so called transverse, longitudinal, polarization and 
interference cross sections, respectively. 
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3.6 Results for eN cNtt 

In fig. 3.5 we show the W = ^/s dependence of the longitudinal (lower line) and 
transverse (upper line) cross sections summing over the final charge of the pion and 
integrating over Jl*. The peak position, the strength and the shape of our results are 
in good agreement with the data of [0]. In general terms our results are very close to 



those obtained in [||] with a similar good agreement with experiment as found there. 

In fig. 3.6, and as an illustration of our discussion about gauge invariance in sect. 
3.3, we show the differences between the results for the inclusive cross sections ut and 
(Tl, obtained: by using Fa and taking (as gauge invariance forces) Ff = F^^^t, = Fa 
(dotted lines); by using F^tttt and taking Ff = Fa = F^^^t, (dash lines); by using Ff and 
taking Fa = F^^.^, = Ff (full lines) and, finally, by using Fa, -F^tttt and Ff (dash-dotted 
lines). As one can see, the differences are negligible, relative to present experimental 
errors. 

In fig. 3.7 we show the results for daj/d^l./ (sin 9l.^/2) in the channel ep — >• enTr"*" 
and compare them with the experimental data of |]62[ . We can see that the agreement 
is reasonably good and so is the case for dap/dQ^/ sin^ 6** shown in fig. 3.8, where the 
data are again from . 



In fig. 3.9 we show the $* dependence of the cross section in the channel ep emi~^ 
for three different kinematics. The data are from and we see again a reasonable 
agreement with experiment. 

Finally in fig. 3.10 we show the $* dependence for another channel, the ep —>■ epn^, 



in order to show a case where the agreement with the data, in this case from ref. [£3|, is 
not as good as in general terms. Given the fact that the angle will be integrated in 
the (e, e') reactions in nuclei, such punctual discrepancies will not matter in our study 
of the nuclear processes. 
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Fig. 3. 5 W = ^/s dependence of the longitudinal (lower line) and transverse (upper line) cross 
sections summing over the final charge of the pion and integrating over fi* . Experimental 
data from [|6lll. 
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Fig. 3. 6 Inclusive cross sections ax and aL obtained: by using Fa and taking Ff = F^j^^^ = Fa 
(dotted lines); by using F^^^^^ and taking Ff = Fa = -F^tttt (dash lines); by using Ff and 
taking Fa = F^^^-,^ = Ff (full lines) and, finally, by using Fa, F^t^^^ and Ff (dash-dotted 
lines). Experimental data from []6T| . 
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Fig.3.7 Calculation of da i /dni/ {sin 61^/2) in the ep — > envr^ channel . Experimental data 



from [32|. 
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Fig. 3. 8 Calculation of dapjdVl*:^! sin^ 0* in the enir^ channel. Experimental data from 
32i. 
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Fig. 3. 10 Calculation of the <I>* dependence of the cross section in the ep — > enir^ channel. 
Experimental data from |63]. 

In the next sections we shall use this model to evaluate the pion production con- 
tribution to the (e, e') cross section, as well as the exchange currents which contribute 
to the 2N emission channel. 



4 The (e, e') reaction in nuclei. 
4.1 Formalism 

We want to use a covariant many body formalism to evaluate the (e, e') cross section. 
For this purpose we evaluate the electron self-energy for an electron moving in infinite 
nuclear matter. Diagrammatically this is depicted in fig. 4.1 
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n 

Fig.4.1 Diagrammatic representation of the electron self-energy in nuclear matter. 

The electron disappears from the elastic flux, by inducing Iplh, 2p2h... excitations 
or creating pions, etc., at a rate given by 

r{k) = -2^/mS. (20) 

where ImE is the imaginary part of the electron self-energy. This latter magnitude 
can be readily evaluated from the diagram of fig. 4.1 and we find: 

^^^^^ = J (^"-(^)7. fe.2 _ ^2 ^ ie ^-"-^^^ (7W ^ ^ 

where H^'^ is the virtual photon self-energy. Eq.(21) displays explicitly the electron 
propagator (fraction after 7^) and the photon propagator (5^ + ie)~^ which appears 
twice. By averaging over the spin of the electron, r, we find 

'^^''^ " 2mJ (27r)^ {k" -ml + ie) ^ ' 

and since we are interested in the imaginary part of E we can obtain it by following 
the prescription of the Cutkosky's rules. In this case we cut with a straight horizontal 
line the intermediate e' state and those implied by the photon polarization (shaded 
region). Those states are then placed on shell by taking the imaginary part of the 
propagator, self-energy, etc. Technically the rules to obtain /mS reduce to making the 
substitutions: 

E(A;) ^ 2iImE{k)Q{k") 

E(k') ^ 2i/mS(A;')e(A;'°) (23) 
n^^^(g) ^ 2ilmn'"'{q)e{q^) 

where 

= w2 (24) 
k — m^ + ie 

and © is the Heaviside, or step, function. By proceeding according to these rules we 
obtain 
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The relationship of JmS to the (e, e') cross section is easy: TdtdS provides a prob- 
abihty times a differential of area, which is a contribution to a cross section. Hence we 
find 

2/77 '2,771 

da = T{k)dtdS = -—ImEdl dS = -^/mSrfV (26) 
and hence the nuclear cross section is given by 

/ 2,171 
d\—ImE{k,p{f)) (27) 
\k I 

where we have substituted S as a function of the nuclear density at each point of the 
nucleus and integrate over the whole nuclear volume. Eq. (27) assumes the local den- 
sity approximation, which, as shown in [^, is an excellent approximation for volume 
processes like here, hence we are neglecting the electron screening and using implicitly 
plane waves for the electrons (corrections to account for the small distortion are usually 
done in the experimental analysis of the data, see 
Coming back to eq. (25) we find then 



d^a 


a \k'\ 1 


dVt'JE'^ 


\k\ (2vr)2 



j d\ (/mn^;^ (2^ 



which gives us the (e, e') differential cross section in terms of the imaginary part of the 
photon self-energy. 

If one compares eq. (28) with the general expression for the inclusive (e, e') cross 



L^'W^, (29) 



section |7T|, |72| (see also eq. (11)) 

d^a c? \k'\ 



^4 I ri 



we find 



W^''' = ^ / d^r- {ImW + ImW^ (30) 

ne J 2 

Once again, by choosing q in the z direction and using gauge invariance one can 
write the cross section in terms of the longitudinal and transverse structure functions 
Wl, Wt as 

-^^"^ ( + (31) 



dVt'JE'^ \dn 
where 



Mott 
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da 

m 



Mott 



a^cosl{e/2) 
AEls\nl{e/2) 



(32) 



and 



(33) 



Wt 



Hence using eq. (30) we can write Wl and Wt in terms of the photon self-energy 



where we see that we only need the components II™ and 11^^. 

4.2 The virtual photon self-energy in pion production 

We must construct a self-energy diagram for the photon which contains pion production 

in the intermediate states. This is readily accomplished by taking any generic diagram 
of the 7*A^ — > ttN amplitude of fig. 3.2 and folding it with itself. One gets then the 
diagram of fig. 4.2 where the circle stands for any of the 6 terms of the elementary 
model for 7*A^ — > nN. The hues going up and down in fig. 4.2. follow the standard 
many body nomenclature and stand for particle and hole states respectively. 



Fig.4.2 Photon self-energy obtained by folding the 7*iV — ttN amplitude. 

The photon self-energy corresponding to this diagram (actually 36 diagrams) is 
readily evaluated and gives 



as 




(34) 
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d^k (fp nN{p)[l - nN'{p + q - k)] 

{27i'f J [27^)^ - k"" + E{p) - E{p + q - k) + le^ ^35^ 



xD^{k)lTr^P'''{T^'T^^] 



NN' 



where is the amphtude for 7*A^ ttN. The indices A^, A^' in eq. (35) stand for the 
hole and particle nucleon states respectively and rijqij)) is the occupation number in 



the Fermi local sea. E{p) is the energy of the nucleon yp^ + and is the pion 
propagator 

D.{k) = -2 V— - (36) 

k — + It 

A further simplification can be done by evaluating the amplitudes at an average 
Fermi momentum. Explicit integration over p and also this approximation were done 
is and the approximation was found to be rather good. We take < p >= \J\kp 
with kp the local Fermi momentum (37r^p(r)/2)^/'^ and a direction orthogonal to that 
of the virtual photon. The errors induced by this approximation are smaller than 5%. 
Then we can use the Lindhard function Un,n' defined as 

U ia k)-2f "^'^ nr{p)[l-ns{p + q-k)] 

"''^ ' J {27^fq'^-k^ + E{p)~E{p + q-k)+ie ^ ' 

where the indices, r, s correspond to protons or neutrons. 

For the evaluation of the imaginary part we need an extra Cutkosky rule 

U{p) 2iImU{p)Q{p^) 

which is the general rule considering that the Lindhard function plays the role of a ph 
propagator. 

Hence we apply the Cutkosky rules of eq. (23) and the former one and we find 

ImW'^= f^ImU^^iq-k)-^e{q'-u{k))x 

J (2vr) 2u:{k) ^^^^ 



NN' 



kO=uj{k) 



Since there are analitycal expressions for IttiUnn' (see Appendix B of [§^|), the ap- 
proximation done saves us three integrals and a considerable amount of computational 
time. 

There is an interesting test to eq.(38). Indeed, in the limit of small densities the 
Pauli blocking factor 1—n becomes 1 and ImUNN'iq) — "~^Pa^^(Q'°~ 2M^)- Substituting 
this into eq. (38) and (29), (30) one easily obtains that aeA = <^epZ + cfenN, the strict 
impulse approximation. By performing the integral in eq. (38) one accounts for Pauli 
blocking and in an approximate way for Fermi motion. Later on we shall introduce 
other corrections due to medium polarization. 
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4.3 The A excitation term 

One of the terms implicit in eq. (38) is the one where one picks up the A excitation 
term both in and is T^". This term is depicted diagrammatically in fig. 4.3(a) 
and, hke in pion-nuclear and photo-nuclear reactions at intermediate energies, plays a 
major role in this reaction. 

In order to evaluate this piece one can go back to eq. (35) and perform the d'^k 
integration to factorize the A width and on the other hand one will also have the 
modulus squared of the A propagator present in the A term of eq. (5). This, however, 
can be obtained more economically by reinterpreting the diagram 4.3 (a) as a Ah 
excitation with a A width. We can also divert a little from the general formulation, 
and in order to gain some extra accuracy we can implement Lorentz covariancc exactly 
simply boosting the tensor U^^ from a frame where the A is at rest {q + p = PA = 0), 
where the amplitude of eq. (5) would be (by construction) equivalent to the relativistic 
amplitude. 






Fig.4.3 Diagrammatic representation of the Ah photonuclear excitation piece. 
Hence we get 



12^7 



xTr ((5t X q^mnS X q^mY) ^ 



■X 



(39) 



/r»S'\(p + r/) - -{p + q) 
X 2 



V^-MA + i^^-EA(s) 



The Lorentz matrix A is such that A-y^g^ = q^^ and s — {p + q^. The coefficients 
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Cij, i,j = 1, 2 stands for proton or neutron, account for isospin and are given by 



= o + + (40) 



'2\2 2, 2 

Eq. (39), however, seems to neglect the Pauh blocking factor 1 — n of eq. (35). This 
factor, however, is taken into account implicitly and leads to the Pauli blocked width 



r, which is evaluated in ||6^. Furthermore, in a nuclear medium the A is renormalized 
and acquires a self-energy Sa, which is also accounted for in eq.(39). The results of 
|51| for Ea are used in the calculation. This self-energy accounts for the diagrams 



depicted in fig. 4.3, where the double dashed line stands for the effective spin -isospin 
interaction, while the serrated line accounts for the induced interaction. The effective 
spin-isospin interaction is originated by a pion exchange in the presence of short range 
correlations and includes p-exchange as well. It is obtained by substituting 

mjDM) miVM) + {^ij - mj)yt{q) (4i) 

and expressions for VJ, Vt are found in ( {f-n-NN/fn-KyVi,t here is equivalent to Vi^t of 



|5T|). The induced interaction accounts for the series of diagrams depicted in fig. 4.4. 



Fig. 4. 4 Feynman diagrams included in the evaluation of the A self-energy. 

There is an RPA sum through ph and Ah excitation and is readily obtained as 

V. , - an 



l-Uiq)Vi{q) 
{Sij-qiqj) 



NN 



Vtiq) 



1 - U{q)Vt{q) 



(42) 

"I 



NN 



where now U{q) = UN{q) + U/\{q) is the Lindhard function for ph + A/i excitations 
including forward going and backward going bubbles ||^ in contrast to U which only 
contains the forward going bubble of a ph excitation (the only one which contributes 
to IuiUn for g° > 0). in addition incorporates a factor two of isospin with respect 
to U, such that ImUN = 2ImU for symmetric nuclear matter. However, all the work 



which goes into the evaluation of Sa is done in ref. fSj], where a useful analytical 
parameterization of the numerical results is given that we use here. The imaginary 
part is parametrized as 
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(43) 



where the different coefficients are given in |]^ as a function of the energy. 

The separation of terms in eq. (43) is useful because the term Cq comes from the 
diagrams (c) and (d) of fig. 4.3 when the lines cut by the dotted line are placed on 
shell, and hence the term is related to the (7*,vr) channel, while Ca2i Caz come from 
the diagrams (b) and (e) and are related to two and three body absorption. Hence, 
the separation in this formula allows us to separate the final cross section into different 
channels. 




Fig. 4. 5 Irreducible pieces in the A/i channel from the A/i interaction. 

It is also easy to realize that the RPA sum of Ah excitations, shown in fig. 4.5, can 
be taken into account by substituting ReL/^ by |^| 




(44) 




Fig. 4. 6 Diagrammatic representation of the inclusion of a ph excitation between A/i excita- 
tions. 

and furthermore, if we wish to include some ph excitation in between, see fig. 4.6, 
(which is actually not relevant numerically), this is done easily by substituting i^eS^ 
by 



If 11 
9 I 



f2 

JttNN 



(45) 



UNVt 
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4.4 Results for the A contribution 



In fig. 4.7 we show the results coming from the A term discussed in the former section. 
The experimental data are coming from ||^. We have separated the contribution from 
the different channels. Besides the upper solid line which stands for the total contri- 
bution, looking from up to down at about uo = 350 MeV the next line corresponds to 
pion production, the following one is two nucleon absorption and the lowest one three 
body absorption. We can see that most of the experimental strength in the A region 
is provided by this A excitation term, but there is still some strength missing. In fig. 
4.8 we show the contribution of the delta piece for the "^^^ Pb nucleus. The data are 
now from IffSl. The results are similar to those found in ^^C and there is still some 



strength missing. The study of this missing strength will occupy the next sections. 

6 
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Fig. 4. 7 Contribution of the A piece to the (e,e') cross section in ^^C. Experimental data 
from [Hi. See text for different contributions. 
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Fig. 4. 8 Contribution of the A piece to the (e, e') cross section in ^'^'^Pb. Experimental data 
from [73|. 
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4.5 Two body photoabsorption 



Let us go back to the generic diagram of pion electroproduction of fig. 4.2. Let us take 
the pion hne and allow the pion to excite a ph. This leads us to the diagram of fig. 



Fig. 4. 9 Photon self-energy obtained from the one in fig. 4.2 when the pion is allowed to 
excite a ph. 

This is still a generic diagram which actually contains 36 diagrams when in the 
shaded circle we put each one of the terms of the 7*A^ — > nN amplitude of fig. 3.2. 
One must avoid the temptation of factorizing these amplitudes in order to evaluate 
these diagrams since some of them might be symmetric and then have a symmetry 
factor 1/2. This is the case here with one diagram implicit in fig. 4.9, which is the one 
corresponding to the pion pole term in each one of the 7*A^ — > ttN amphtudes. This 
diagram is shown explicitly in fig. 4.10. 



4.9. 





Fig. 4. 10 Pion pole term included in fig. 4.9. 

One can see that the diagrams in fig. 4.9 contribute to /mil according to Cutkosky 
rules when we cut by a horizontal line and the 2p2h are placed on shell. 
The contribution of the diagram of fig. 4.9 is readily done. We obtain 




(46) 



'TT 
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where Ux is the Lindhard function for ph by an object of charge A: this is, twice f/p^„ or 
Up^n for the excitation by a charged pion or Up^p + ?7„^„ for the excitation by a neutral 
pion and k is the pion momentum. The factor F^{k), where is the pion form factor 
appears because now the pions are off shelL Recall that we also take all form factors 
equal in order to preserve gauge invariance (eq. (9)). The factor Sa is the symmetry 
factor, unity for all diagrams and 1/2 for the symmetric one of fig. (4.10). 

We can again simplify the expression by taking an average nucleon momentum of 
the Fermi sea to evaluate the matrix elements of T^T'^^ . This allows us to factorize 
the Lindhard function and we get 



NN' ^ ' 1" 

PUx{k)\Tr^P^%T^^T^'')NN'S^F^{k) 
By applying Cutkosky rules we find 



(47) 



NN' 



(Tk 



e(g° - k^)ImUNN'{q - k)e{k°)x 



ImUx{k)Dl{k)i^pF^{k)S, 



(4J 



INN' 



The cut which places the two ph on shell in the diagrams of fig. 4.9 is not the only 
possible one. In fig. 4.11 we show a different cut (dotted line) which places one ph and 
the pion on shell. 




Fig. 4. 11 Same as fig. 4.9 and showing the cut which places one ph and the pion on shell. 

This contribution is taken into account in the A excitation term by means of the 
term Cq. As done for real photons in [Q, we neglect this contribution in the other 
terms, because at low energies where the background pieces are important, the (7*, vr) 
channel is small and at high energies where the (7*,vr) contribution is important, this 
channel is dominated by the A excitation and there this correction is taken into account. 
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We have also considered two body diagrams where each photon couples to different 
bubbles: As found is only one of them is relevant, the one in fig. (4.12), which 
involves the KR term alone and which we take into account. 




Fig. 4. 12 Feynman diagram related to the KR term of fig. 4.9 with outgoing photon from 
the second nucleon. 

Following the same rules as above, this term is readily evaluated and gives 

2 



/mil 



00 



(27r) 



{kq){k - q)q. 



N 



(49) 



ImUp.p{q)ImUpp{k - q) + ImUpn{q)ImUnp{k - g) + 



+ ImUnp{q)ImUpn{k - q)] 0(A;O)0(A;O - q^) 



/ ^Mk)Do{k - q)\k\htn^eF^^^{k)F^^^{k - q)x 
\ (t^-) ImUpp{q)ImUpp{k - q) + 
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2M 



V2 



2M ) ^ 



ImUpn{q)ImUnp{k - q) 
ImUnp{q)ImUpn{k - g) ^ x 



(50) 



e(A;0)e(fcO-gO) 

The contribution of this term is roughly 1/2 of the KR x KR term in the generic 
diagram of fig. 4.9. 

4.6 Contributions tied to the (7*,27r) channel 

The — > ttttN reaction has been the subject of recent detailed experimental analyses 
Iff^ , 75 1 and also of recent theoretical analyses, some of them spanning a large energy 
range []76| , [77| and others concentrating only very close to threshold in order to test 
predictions of chiral perturbation theory ||78|, [79| . 

The model in ref. ||76|] for the 7p -k^-k'p uses 67 Feynman diagrams, while ref. 



|77| , where the model is extended to the other isospin channels uses only 20 diagrams 
which are necessary below = 800 MeV, where the new data have been measured. 

Although the model is rather elaborate and contains many terms, one can see that 
the gross features of the reaction can be obtained with the two terms of fig. 4.13, 
accounting for about 80% of the cross section. 




Fig. 4. 13 Relevant Feynman diagrams that enter in the evaluation of the 7*A^ N2tt cross 
section. 

Since here we are only concerned about corrections to the more important terms 
which we have discussed above, it is sensible to just take these two diagrams. The 
diagram in fig. (a) is the ANir'y Kroll Ruderman (KR) term, while the one in fig. (b) 
is the pion pole term. In both terms the A is excited. The KR term, which appears 
from the AA^vr vertex by minimal coupling, is given by. 



7T^ 




(51) 
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corresponding to tt^A production in the '-f*N ^ Att vertex. 

By following the same steps as before we obtain the many body diagrams of fig. 
(4.14), where the dashed circle indicates any of the two terms of fig. 4.13. 




Fig. 4. 14 Photon self-energy diagrams obtained by folding the terms of fig. 4.13. Diagram 
(b) is obtained when the pion is allowed to produce a ph excitation. 



Furthermore, as discussed in ref. in the diagrams of fig. 4.14 (b) we keep only 



the term with the KR x KR in the vertices. This is done since the pion in flight term, 
can be considered as a two step process of a 7*A^ —>■ ttN with tt, a real pion, followed 
by the ttN — > A excitation. The two step processes redistribute strength but do not 
change the cross section and hence are not included in our approach. 

Since the A in the Ah excitation in fig. 4.14 is also renormalized, we are accounting 
for the physical channels depicted in fig. 4.15 when placing on shell the states cut by 
the dotted line. 




(1) ^ (2) (3) -J" (4) 

Fig.4.15 Detail of fig. 4.14 indicating the physical channels associated to the cuts. 



As one can see there, (1) accounts for lpl/i27r excitation, (2) and (3) for 2p2/il7r 
excitation and (4) for 3p3h excitation. 

The evaluation of these pieces follows exactly the same steps as for figs. 4.2 and 
4.9, simply replacing the 7*A^ — > nN by the 7*A^ —>■ ttA amplitudes and one nucleon 
propagator by the A propagator. The contribution of these terms below u) = 350 MeV 
is very small. Their importance increases with the energy and at a; = 450 MeV they 
account for about 1/5 of the cross section, as found for real photons. 
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4.7 Polarization (RPA) effects 



In the diagrams of fig. 4.9 we can consider the ph as just the first order of a series 
of the RPA excitations through ph and Ah excitations. If one replaces the ph by the 
RPA series, one is led to the terms implicit in fig. 4.16. A similar series would appear 
for the case of the (7,7r) process depicted in fig. 4.2. 

In practical terms this is done in a simple way by having a bookkeeping of both the 
spin longitudinal and spin transverse parts and replacing 



ImU, 



ImU, 



N 



N 



ImU, 



N 



\l-Ux{q)Vi\^ 

where a, h measure the strength of the longitudinal and transverse parts. 



(52) 



For the transverse part of the photon self-energy , 11^^ the procedure to follow is 



identical to the one explained in section 9 of and we refer the reader to this paper 
(see also ||5^). The only novelty here is 11'^'', but this component is of spin longitudinal 

,6 = 0. In the A term 
hence, no 



character and is renormalized by means of eq.( ^) with a = 1, 6 = 0. 
the polarization effects are already included in the self-energy of ref. 
further corrections are needed. 
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5 Short range correlations 

So far the calculations have been done using implicitly plane waves for the nucleon 
states. Short range nuclear correlations modify the two nucleon relative wave function 
and this has a repercusion in some of the matrix elements which we have calculated. 
This is particularly true in those matrix elements which involve a p-wave coupling in 
the vertex for each of the two nucleons, because the pion exchange generates a 5{r) 
function, which is rended inoperative in the presence of short range correlations. This 
is not exactly the case if finite size effects by means of form factors are taken into 
account, but the need to implement the effects of the short range correlations remains. 
The correlations can also introduce spin transverse components in the p-wave-p-wave 
terms which were originally of the spin longitudinal nature |^5|. Hence, at the same 



time that one introduces the effect of correlations, one takes advantage of this and 
introduces the p meson exchange in this case. In this way, we generated Vi and Vt of 
eq.(41). 

The method to introduce the effects of correlations is to substitute a two nucleon 
amplitude V{q) by 



{2ny 

where ^l{p) is the Fourier transform of a nuclear correlation function. 

Once again the techniques to make these corrections can be seen in (appendix 
D, see also [H]), and we do not repeat them here. 
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6 Detailed study of the quasielastic peak 



6.1 Formalism 

So far we have studied the (7*,7r) process and the 7* absorption by a pair or trio of 
particles. This was done keeping a parallehsm to the real photon case. However, unlike 
the case with real photons, a virtual photon can be absorbed by one nucleon leading to 
the quasielastic peak of the response function. We have left this problem till the end in 
order to introduce the concepts of many body which proved relevant in the scattering 
of real photons with nuclei and in the equivalent channels of virtual photons studied 
before. 

Then we can use the same concepts and ideas here in order to introduce the appro- 
priate many body corrections to the quasielastic peak. 

Thus we begin by evaluating U^^^ for the Iph excitation driven by the virtual photon, 
as depicted in fig. 6.1. 




Fig. 6.1 Photon self-energy diagram for the Iph excitation driven by the virtual photon. 



The photon self-energy associated to this diagram is given by 
"7(2^ p"-E{p)+ze'' 



(53) 

X ^(l-n,■(p + g)) T^ryi^v^-) 

p" + - E{p + q)-ie ^ ' 

where represents the ^NN vertex and E{p) is the nucleon kinetic energy. The 
vertex is given by 

Once again the application of Cutkosky rules leads to 

/mn™ = -/^ ni{p){l-nj{p + q))x 

■' ^2^) (55) 

x5{f + g° - E{p + q) )Tr( WtO) 



36 



ni{p){l -nj{p + q))x 

x5{p° + g° - E{p + q) )Tr(K^0^) 



(56) 



and if we average Tr{V'^V'^'^) over the nucleon momentum in the Fermi sea we can 
write 



/mn°° = -ImU{q,p){Tr{V°V^°)) 



(57) 



ImW = ^/mf/(g,p)(rr(r^0^)) 



(5^ 



The average over the Fermi momentum can be done keeping terms up to q^ /M^ 



N 



and we find in terms of 



A 



00 



1 - 



4M^ 



G|(g) 



GiM) 



\i:^v' + cff-\fGlAq)y 



(59) 



Ml 



1 - 



N 



- \q'GlM) \ (60) 



4M 



N 



where p^ = + and Ge,Gm are the Sachs form factors pOl, ITT . 
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6.2 Spectral function description and final state interaction 

One of the corrections to the bare ph excitation studied above is the one induced by 
final state interaction, as we indicated in section 2, which in our approach can be taken 
into account by dressing up the nucleon propagator of the particle state in the ph 
excitation, as depicted in fig. 6.2 (there the dashed line would account for the whole 
NN interaction not just pion exchange). However, some caution must be exerted when 
talking about this diagram. 
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p+q 



Fig. 6. 2 Photon self-energy diagram obtained from fig. 6.1 by dressing up the nucleon pro- 
pagator of the particle state in the ph excitation. 

In the first place, this is one of the terms implicit in the generic diagram of fig. 4.9 
when the nucleon pole term is taken in each of the 7*A^ — > A^A^ amplitudes. This term 
poses no problem for real photons and leads to a small fraction of the two nucleon 
absorption. However, for virtual photons this diagram is divergent. The reason is that 
when placing the 2p2h excitation on shell through Cutkosky rules, we still have the 
square of the nucleon propagator with momentum p + q in the figure. This propagator 
can be placed on shell for virtual photons (not for real photons) and we get a divergence. 




Fig. 6. 3 Insertion of the nucleon self-energy on the nucleon line of the particle state. 

The divergence is physical in the sense that its meaning is the probability per unit 
time of absorbing a virtual photon by one nucleon times the probability of collision of 
the final nucleon with other nucleons in the infinite Fermi sea in the lifetime of this 



nucleon |£6|. Since this nucleon is real, its lifetime is infinite and thus the probability 
infinite. The problem is physically solved by recalling that the nucleon in the 
Fermi sea has a self-energy with an imaginary part which gives it a finite lifetime 
(for collisions). This is also immediately taken into account technically by iterating 
the nucleon self-energy insertion of fig. 6.3 on the nucleon line, following the Dyson 
equation, hence substituting the particle nucleon propagator by a renormalized nucleon 
propagator including the nucleon self-energy in the medium. 
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G(p°,p) = ^ (61) 

where J2{p^yP) is the nucleon self-energy. Alternatively one can use the spectral func- 



tion representation ^7 



Gip'p, = /' f"<"-'^ + r 4^^d. (62) 

^ ' J~oo -uo-ie Jf, -uo + it ^ ' 



where Sh, Sp are the hole and particle spectral functions related to S by means of |^ 
*Ifa;>/z , Spiu;,p) = -^ImG{u,p) = -l ^^f^Y^ 



*Ifc.</i , S,{uj,p) = hmG{uj,p) = y^^i^^P) 



(63) 



and /i is the chemical potential and 



2 N 2 



B = {ImTi{uj,p))'^ 

By means of eq.(62) we can write the ph propagator or new Lindhard function 
incorporating the effects of the nucleon self-energy in the medium and we have for 
ImU 



ImU{q) = -7j- / dpp^ dx dujSh{uJ,p)x 

JO J-l Ju-q" 



(64) 



xSp{q^ + uj,\/p +q +2pqx) 



We use the spectral functions calculated in [^] , but since the imaginary part of the 
nucleon self-energy for the hole states is much smaller than that of the particle states 
under consideration we make the approximation of setting to zero /mE for the hole 



states. This was found to be a good approximation in |45]. Thus, we take 



S,{u;,p) = 6{u; - E{p))Q{fi - E{p)) (65) 

where E{p) is the energy associated to a momentum p obtained selfconsistently by 
means of the equation 

E{p) = ^ + RenE{p).P) (66) 
The chemical potential was then taken as 
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2M, 



AT 



+ Re'E{ii, kp) 



where kp is the Fermi momentum. 

It must be stressed that it is important to keep the real part of E in the hole states 
when renormalizing the particle states because there are pieces in the nucleon self- 
energy largely independent of the momentum and which cancel in the ph propagator, 
where the two sclfcncrgics subtract. 

The effect of the use of the spectral function, accounting for FSI is a quenching of 
the quasielastic peak and a spreading of the strength at higher energy as can be seen 
in fig. 6.4. 



c; 

t3 




120 140 
u (MeV) 



220 



Fig. 6. 4 Effect of the use of the spectral function in the evaluation of the Lindhard function. 
The uncorrclatcd Fermi sea results are obtained from eqs.(57),(58). Those with the medium 
spectral function, with the same equations substituting the bare Lindhard function U by the 
medium modified one of eq.(64). 



6.3 Polarization (RPA) effects in the quasielastic contribution 

We take now into account polarization effects in the Iplh excitation, substituting it by 
an RPA response as shown diagrammatically in fig. 6.5. 
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Fig. 6. 5 Diagrammatic representation of the polarization effects in the Iph excitation. 
For that purpose we use an effective interaction of the Landau-Migdal type 

V^(fi, fa) = Co6{fi - fa) {/o(p) + /o(p)n^2 + 



and we take the parameterization for the coefficients from ref. 



(67) 



fMr)) = ^ft^ 



P(0)' 



p(r) 



(ex) 



Ain) 
JO 

Aex) 

Jo 

/{in) 



0.07 

-2.15 
-- 0.33 



P(0)_ 
).45 

Co = 380 MeVfm^ 



(68) 



/?^") = 0.45 



fo 

gt^ = gt^ = go 
g't^ = g't^ = g'o = 0-725 



(69) 



0.575 



For consistency, in the S = 1 = T channel {aarr operator) we have continued to 



use the interaction used in 32 which has been used for the renormalization of the 



pionic and pion related channels studied in the former sections. There is only a minor 
difference of about 4% in q'q between the two parametrizations. 
Recalling that we had 



(70) 



let us take the nonrelativistic reduction of A'^'^ in order to see the effects of the RPA 
renormalizations 



N 



+ (p n 



(71) 
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Given the spin-isospin structure, the electric and magnetic components will be 
renormalized in the following way: 

a) Interaction aarf: is the one we used to renormalize the pionic related channels 
in former sections. It affects only the magnetic components. If we write 



mag. 



Xt 2 — 



(72) 



+ [neutrons] 2^'^'^ 



it is easy to see that the magnetic part of II*-' becomes 

JttI 



mag. 



\UNA'iag.{q)e' + lf^^ 



4M^ mi ^ 



Vt{q) 



Vt{q)U{q) 



(73) 



'5'^ -q'q')GlM) {f^pUp-f^nUr. 



where U = 11^ + Ua- 
b) Interaction rr: 

This interaction selects the non magnetic components of V^. Thus and the 
convective terms of A^^ (term with 2p + g in eq.( ^)) are renormalized. 

However, given the smallness of the convective terms (about 10% contribution to 
the transverse response) we shall not consider their renormalization. 

Thus we consider only the modification to A^^ from this source. Since A^^ is given 

by 



^00 



+ 



J2XrF?{q)X'Xr'F?{q)Xr 



Y.XrF^{.q)Xr'Xr'F^{.q)Xr 



the renormalized expression for 11°° will be 



— + 



2 



(74) 



noo ^ g2 1 ^ppy^^ ^ ^pny^^ ^ Cofl^iFfU, F^U, 



riTT N2 ■ 



1 - CofoUNiq) j ^^^^ 

where in the denominator we do not have now 11^ since the operator fr cannot excite 
A components. 

c) Interaction aa: 

Here again, like in case a), only the magnetic components are modified. We find 
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X (g 25^- _ q^q^)G\M) + /i„f/„) ^^^^ 

+ e//ect of (a) 

The correction from the RPA sum is taken into account by means of the second 
term of the right hand side of eq.( |7^) to which we should add the same term in eq.( ^3|) 
coming from the renormahzation with the aaff operator. 

d) Scalar interaction: 

This one affects A^^ and we find 

= \{F!fu, + {F-fu^ + ^Mp){f!u +F^U^f \ ^ ^^^^^ ^^^^ 

[ t - cojoUNyq) ) 

We show in figs. 6.6, 6.7, 6.8 and 6.9 the effects on Rl and Rt {Rl = —-^-t'Wl 

and Rt = 2Wt) of the different polarization terms. The solid line corresponds to 
the calculation including these effects and the dashed line, to the calculation without 
polarization effects (we are using spectral functions in the calculation of the Lindhard 
function): 

•1 I I I I I I I I I I I I ^ I I I ^ I I I ^ I I I ^ I I I ^ I I I I I I I 



Scalar Interaction 
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Fig.6.6 Polarization (RPA) effect (solid line) in the evaluation of Rl: scalar interaction. 
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Fig.6.7 Polarization (RPA) effect (solid line) in the evaluation of Rl- tt interaction. 
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Fig.6.8 Polarization (RPA) effect (solid line) in the evaluation of Rt- So interaction. 
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Fig.6.9 Polarization (RPA) effect (solid line) in the evaluation of Rt'- Butt interaction. 
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One possible source of renormalization not yet considered is the one shown in fig. 
6.10, where a ph is attached to one photon and a Ah to the other one, plus any other 
possible ph or Ah excitations in between. 




Fig. 6. 10 Diagrammatic representation of a possible source of renormalization: a ph is at- 
tached to one photon and a to the other one, plus any other possible ph or Ah excitations 
in between. 



This contribution affects the transverse part and both the quasielastic as well as 
the A peak. However, given the small interference between ph and Ah excitations, the 
contribution of these terms is not significant. We give, however, the expression here 
for completeness 

n«i - e f 4f A 2 

-^U; m^l r Y 



The effect of the polarization is moderate, but relevant when aiming at a precise 
description of the process. 



We show in figs. 6.11 and 6.12 the effects of the polarization in the longitudinal 

and transverse response functions. The net effect in the cross section is a quenching in 
the quasielastic peak and a spreading of the strength at higher energies. 
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Fig. 6. 11 Polarization (RPA) effect in the evaluation of Rl. 
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Fig.6.12 Polarization (RPA) effect in the evaluation of Rt- 

6.4 Further considerations 

Some other terms appearing in the generic diagram of fig. 4.11 require some special 
thought. These are the terms in which one of the vertices contains the nuclcon pole 
term of the 7*A^ — > vrA^ amplitude, while the other one contains all terms but that one. 
This is depicted in fig. 6.13. 
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Fig. 6. 13 Photon self-energy diagrams in which one of the vertices contains the nucleon pole 
term of the 'y*N ^ nN amplitude, while the other one contains all terms but that one. 

Unlike the case of real photons where only the cut (a) exciting 2p2h gives rise to an 
imaginary part, now the cut (b) placing the ph on shell is a source of imaginary part 
which produces strength is the quasielastic peak. The diagrams in fig. (6.13) for the 
cut (b) could be then considered an ordinary ph excitation with a renormalized vertex. 
We have evaluated the two sources of imaginary part in U^'' and find for the cut (a) 



ImU^'^ = i^Y: J A^n,{m^Mp')S / - E{p) - Re^l^,p) x 



5 (pO + - - E{p + q-k)- ReT. + g - fc)) x 



N 



x(l - nj{p + q- k))- 



ImUxjk) 



1 - 



ml 



Vi{k)U{k) 



'l-ni{p + q))x 



x2Re < 



+ - E{p + g ) - E 



N 



KR,PP,NPC,A,AC 



jp + q) 

2Mn 



X 



,p + q 



(79) 

where Mnp is the nucleon pole amplitude of 7*A^ — > vrA^ omitting the nucleon propa- 
gator. On the other hand the contribution of the cut (b) is given by 
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p ) I X 



N ' 



l-g'\^\ U{k) 



^0) 



xe(/)(27r)2(l - ni{p + g))e(/ + q^)x 
x6 (^p'> + q^-E{p + q)-ReJ: + x 

Jj2^''^' \pO + qO-k^-E{p + q-k)+te 
^ (a;^ - - n A;^ - + ze) ^ 

where 11 is the pion self-energy in the nuclear medium 

n = ( (81) 

/ / JnNN 



the subtraction of the free pion propagator appearing in eq.(80) guarantees that in the 
limit p ^ the correction to the jNN vertex vanishes as it should be. 



6.5 Considerations on gauge invariance 

Most of the theoretical models found in the literature on inclusive (e, e') scattering 
from nuclei do not preserve gauge invariance. The requirement of invariance under 
gauge transformations leads to relations between the components (charge and spatial 
current) of the hadronic current which determines the nuclear response. Actually, 
the longitudinal (charge) multipoles are related to the two transverse (spatial current) 
multipoles of the electric type [BOl. Several prescriptions have been used to restore 



gauge invariance [g^ , ||8T|] . However, in a recent work [^] the arbitrariness of the most 
common prescriptions is discussed in detail with the conclusion that the standard 
procedures to impose gauge invariance in calculations, based on models which do not 
verify it, are misleading and for very low nuclear excitation energies do not ensure at 
all that a better, or a more reasonable, description of the data will be obtained. 

Our model, as we shall see, is gauge invariant at the lowest order (impulse approxi- 
mation) of the density expansion and this symmetry is only partially broken when some 
non-leading density corrections are included. In what follows, we will study the con- 
sequences of the partial breaking of the gauge symmetry for the kinematics studied in 
this paper (from the quasielastic peak to the A excitation region), which involves larger 
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nuclear excitation energies than those studied in [|^. We will also study the feasibil- 
ity of non-gauge invariant models to disentangle between longitudinal and transverse 
channels. 

The unpolarized cross section for inclusive (e, e') scattering from nuclei is given by 
(eq.(29)): 

(82) 



Lorentz, space-inversion, time-reversal and gauge invariance constraint the form of 
the hadronic tensor W^y, which determines the nuclear response. Indeed, the most 
general expression for this tensor assuming the latter symmetries is given by p3[: 



with q and P the virtual photon and initial hadronic system four-momenta respectively 
and M\ = P^. The structure functions 2 are unknown scalar functions of the virtual 
photon variables which determine the nuclear response. 

Using the expression for the hadronic tensor of eq. (|83| ) and taking q in the z 



direction the cross section of eq. (|82D in the lab system becomes (eq. (31)): 
(fa ( d(j\ ( q^ \ \ WT{uJ,\q\) 



WL{uJ,\q\)+ "^'^ (84) 



where the structure functions Wl and Wt are given in terms of Wi and W2 by: 



I -. 1 2 2 2 2 

Wl = -W^i - ^W^2 = -^W^'' = = — ^W^' (85) 



q^ uj^ \q\^ uj\q\ 

Wt = Wi = W^^ = W^y (86) 

Let us now suppose that the hadronic model does not preserve gauge invariance. 
In these circumstances the hadronic tensor (we will call it W^" to differentiate it from 



the one defined in eq. (|83D ) is not conserved (ie, g^W^'^ 7^ 0, q^W^" 7^ ), and is now 



given in terms of four independent functions: 



g2 



pf'p 



A 



+ Ws'-^ + W,—^ (87) 
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Because of the loss of gauge invariance now , and W^^ are no longer related 
and become independent. Thus, we have now: 











Wt = 


Wl = = 


Ws = 






1^1 



(89) 
(90) 



W^ = >V°° + W ~{~ + Jf\j ^9^) 

If gauge invariance is restored and therefore the hadronic tensor is conserved, the 
response functions W3 and W4 vanish and Wl reduces to any of the expressions of 
eq. (pSj). With this new hadronic tensor the differential cross section is now given by: 

-^^"^ ^ + (92) 




with 



Wi(u;,|g|) = lWL{uJ,\q\)-^W,{uj,\q 




I -*I2 



q' / 

4^W°° + 2^W°^ (93) 

The function W3 does not appear in the expression for the cross section because the 
leptonic tensor is conserved. Note that, one can still factor out the differential cross 
section in the form A + B/e and therefore, despite the breaking of gauge invariance, 
one can still compare the results to the experimental response functions obtained via 
the Rosenbluth plot. Then, the breaking of gauge invariance in the theoretical model 
for the nuclear response leads to a redefinition of the response function which has to 
be compared to the experimental one, in the longitudinal channel, and thus one should 
compute Wl given in eq. (0) instead of Wl of eq. (pSf). In fact, the latter one is not 
well defined and there is an arbitrariness in its definition because the 00, Oz and zz 
components of VV"^ are no longer related. 

Traditionally, the longitudinal response function is calculated from the "charge- 
charge" component of the hadronic tensor (— pVV''°). In general, the response 
function calculated in this way will differ from that calculated by means of eq. (0). 
Now we will examine the difference between both approaches as a function of the energy 
and momentum transferred to the nucleus. In order to do that, we define 



5 

.2 

2 



W'' = ^W°° + AW (94) 
>V°" = ^W°° + AW°" (95) 
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where AW^^ and A>V°^ account for the breaking of the gauge symmetry. By construc- 
tion one expects: 



SiT^ (96) 



WOO ^|^|2 

where 5i and 62 would be 1 in the case that AW^ = and AW"^ = >V°''. We 
certainly expect that in our case 61 and 62 are smaller than 1 because at order p in the 
density expansion we exactly fulfill gauge invariance (see comments below). Taking a 
conservative point of view one sees that the ratios in eqs. (96), (97) are at most of 
order 1. Using the definitions of eq. (p4D, we can now write 



The size of the corrections to the "charge-charge" prescription (the term proportional 
to 1 in the formula) traditionally used in the literature depends on the kinematics under 
study. Here we will pay a special attention to three different regions (keeping always 
the momentum transferred to the nucleus smaller than about 500 MeV) : A-resonance 
and quasielastic peaks and the region (dip) between both peaks. 



Quasielastic peak : In this region we have uj ^ |g P/2M7V and then the coefficients 

AW°^ 

yyuir ana yyoo 



of the ratios -^ro- and ^^0^ in eq. ( p6D turn out to be 



^ ^ ^ (99) 



2^ ^ ^ (100) 



g4 M, 



N 



Thus, taking also into account the estimates of eq. (p6D, one finds corrections to 
the "charge-charge" prescription of the order of 62\q\'^ /2M'^ which at most could 
be of the order of 5 — 10% for the momenta and energies transfers studied in this 
paper, assuming 62 ~ 1, which is certainly an overestimate for the reasons pointed 
above. From this discussion, in this region we have decided to use the traditional 
prescription "charge-charge" to compute the longitudinal response function. 

Dip area : In this region the conclusions are similar to those drawn in the previ- 
ous point. However, we would like to point out that gauge symmetry breaking 
corrections to the longitudinal response function are not now as small as before. 

A-resonance peak : In this region uj/\q\ ~ 1 and the situation is radically dif- 
ferent. The corrections due to gauge symmetry breaking are significantly more 
important than in the quasi-free scattering region. For instance, taking the in- 
coming electron energy equal to 620 MeV and the outgoing electron scattering 
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angle equal to 60°, one finds values for the coefficients of the ratios and 
in eq. (|98D of the order of two. Thus the corrections to unity in the bracket 
of eq. (|98|) , are of the order of 2{6i — 62), much larger than in the quasielastic 
peak. We have evaluated Si and 62 in this region and we find Si ~ 0.01 — 0.06, 
62 ^ 0.01 - 0.02 and 2{5i - 82) ^ 0.0 - 0.08. Hence a 10% error in the lon- 
gitudinal response due to the breaking of gauge invariance of our results seems 
realistic in this region. In any case, we should mention that the contribution of 
the longitudinal response to the cross section is very small here and hence the- 
oretical cross sections are largely free of uncertainties due to the small breaking 
of gauge invariance. Since there is no experimental separation of Rl and Rt in 
this region, we do not give these results either. 

We finish this section discussing the origin of the breaking of the gauge invariance 
within our model. Let us consider the diagram of fig. 6.1 whose contribution to the 
virtual photon self-energy in the medium is given in eqs.(53) and (54). One can easily 
check that the imaginary part (when the intermediate nucleons are put on shell) of 
this self-energy is gauge invariant (g^II^'^ oc qfj,V'^ = 0). The first medium correction 
is given by the diagram depicted in fig. 6. 2. When the two ph excitations are put on 
shell, this diagram contributes to the imaginary part of the photon self-energy. This 
new contribution is not gauge invariant because the intermediate nucleon with four 
momentum p + q is not on shell and then the contraction Qf^V^ does not vanish now. 
However, to this level we restore gauge invariance because we consider not only the 
term of fig. 6.2, but also all terms implicit in fig. 4.9. Though the NP amplitude in 
not gauge invariance by itself, the thick dots of fig. 4.9 account for the six amplitudes 
(NP+NPC+KR+PP+DP+DPC) of our model for the eA^ e'Nir reaction. In section 
3, we fixed the different form-factors entering in the amplitudes to end up with a gauge 
invariant model (Eqs.(8-9)). Thus, the leading terms in our density expansion (fig.6.1, 
fig. 4.2 and fig. 4.9) lead to a gauge invariant photon self-energy. 

However, we break again the gauge invariance in section 4.7 when we include the 
polarization corrections to the 36 diagrams of fig. 4.9: we do not renormalize in the 
same way for instance the KR x KR term (which is purely longitudinal and therefore 
gets renormalized only with VJ) than the NP x NP term (which contributes to both 
longitudinal and transverse channels and thus gets renormalized not only with Vi but 
also with Vt). As a consequence, the cancellations in eq.(8) which ensured the gauge 
invariance of the model are altered. 



7 Results 

We have already shown results on the different effects in previous sections. Here we 
will show results with emphasis in comparison with experiment. 

Let us first show results in the quasielastic peak in figs. 7.1, 7.2 we show results for 
Rl and Rt for ^^C and compare them to the data of [17DI. The lower line shows the 



results obtained with the medium spectral function, while the upper one includes also 
the rest of the effects discussed in the former section. 
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Fig. 7.1 Calculation of i?L for C. The lower line in the high energy region corresponds 
to the result obtained with the contribution of the Iplh excitation (fig. 6.3) using the 
medium spectral function of eq.(64). The upper line is the result when one adds the rest of 
contributions: vertex corrections (fig. 6.13), two body absorption diagrams (fig. 4.9), (7*,7r) 
terms (fig. 4.2), (7*,27r) related terms (fig. 4.15), etc. Experimental data from | ]70[ . 
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CO (MeV) 

Fig. 7. 2 Calculation of Rt for ^^C. Same meaning of the lines as in fig. 7.1. Experimental 
data from [^] . 

In fig. 7.3, 7.4, we show results for ^°Ca compared to the data of |^ (fig. (7.3) and 
lower points in fig. (7.4)) and those of [Q (upper points in fig. (7.4)). 
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Fig. 7. 3 Calculation of Rl for ^'^Ca. Same meaning of the lines as in fig. 7.1. Experimental 
data from M. 
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Fig. 7.4 Calculation of Rt for ^'^Ca. Same meaning of the lines as in fig. 7.1. Experimental 
data from (lower points) and [Q] (upper points). 

As one can see, we find a good agreement with the recent reanalysis of 0. 

On the other hand much of the work done here has gone into the evaluation of two 
body mechanisms. In fig. 7.5 we show the results for two body photon absorption 
(solid line) and compare them to pion production (dotted line). Similarly, in fig. 7.6 
we show the contribution of three body photon absorption (solid line) versus the two 
body one (dotted line). We can see that at low energies, the contribution of three body 
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absorption is negligible while at energies around 450 MeV the three body contribution 
becomes sizeable. These results agree qualitatively with the findings of for real 
photons. Let us recall that this classification corresponds to the primary step in the 
collision. The particles produced still undergo secondary collisions in their way out 
of the nucleus. This does not change the inclusive cross section but redistributes the 
strength. The treatment of this FSI and the evaluation of the exclusive channels will 



be treated in a forthcoming paper |£7 
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Fig. 7. 6 Three body photon absorption (solid line) versus the two body one (dotted line) for ^^c. 
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Finally let us see the global results including the quasielastic peak, the dip region 
and the delta region. They can be seen in figs. 7.7, 7.8 and 7.9 for the nuclei of ^^C 
and 20SP6. 

The global agreement is good and the three regions are well reproduced (a bit 
overestimated for ^°®P6). 
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Fig. 7.7 Inclusive (e, e') cross section for ^^C. Ee = 620 MeV and 9e = 60°. The dotted line 

corresponds to the pion production contribution. Experimental data from [p|]. 
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Fig.7.8 Inclusive (e,e') cross section for C. = 680 MeV and 6^ = 36 . Experimental 
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data from [^. 
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Fig.7.9 Inclusive (e,e') cross section for ^°^Pb. Ee = 645 MeV and Oe = 60°. The dotted 
line corresponds to the Iplh excitation contribution. Experimental data from |73| |. 

In fig. 7.7 we also show with a dotted line the results for pion production. 
In fig. 7.9 instead we show with a dotted line the results for the Iplh excitation 
alone. 



8 Conclusions 

We have undertaken the task of constructing a microscopic many body model of the 
(e, e') reaction including all the reaction channels which appear below uj = 500 — 600 
MeV, and which is suited to study the inclusive (e, e') reaction from the quasielastic 
peak up to the A peak, passing through the dip region. Although many studies have 
been devoted to particular energy regions of the spectrum, this is the first work, to our 
knowledge, which ranges this wide energy spectrum. 

Our model has no free parameters. All the input consists of basic couplings of 
photons to nucleons and isobars, and some phenomenological inputs, as correlations, 
which has been tested in former pionic reactions. 

We include explicitly the IN knockout channel, the virtual photon absorption by 
pairs or trios of particles, the pion production plus exchange currents mechanisms tied 
to the {'-/*, 2tt) channel and which contribute to (7*, A^A^vr) or {'-f*,NNN) channels. 

We include effects which have been found important in earlier works, like polariza- 
tion, renormalization of A properties in a nuclear medium, FSI effects through the use 
of spectral functions and meson exchange currents. 

The meson exchange currents are generated in a systematic way from a model for 
the elementary pion electroproduction on the nucleon, which reproduces accurately the 
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experimental data. 

We have payed some attention to the question of gauge invariance, showing that it 
is preserved in our approach in leading order of the density expansion. We also show 
that the appropriate prescription to evaluate the longitudinal response is from the W'^^ 
component of the hadronic tensor, which minimizes the breaking of gauge invariance 
at higher orders in p. 

We evaluate cross sections in the energy range from the quasielastic peak to the A 
peak and find good agreement with experimental data. The three traditional regions: 
quasielastic peak, dip region and delta peak, are well reproduced in our scheme. 

We also separate the longitudinal and transverse response functions in the quasielas- 
tic peak and find good agreement with the latest results of the analysis of Jourdan from 
the world set of data. 

We have used the technique of the local density approximation, which has been 
shown before to be particularly suited to deal with inclusive cross sections and which 
makes unnecessary the use of sophisticated finite nuclei wave functions. 

Finally, the method used here allows the separation of the contribution of different 
channels to the inclusive cross section. This information is the seed to produce exclusive 
cross sections like (e, e'N), (e, e'NN), (e, e'vr), (e, c'-kN) etc. However, this still requires 
to follow the fate of all the particles produced from their production point in the 
nucleus, which is usually done using Monte Carlo simulation techniques, and this will 
be the subject of some future work. 

We would like to acknowledge useful discussions with R.C Carrasco, C. Garcfa- 
Recio and A. Lallena. This paper is partially supported by CICYT contract no. AEN 
96-1719. One of us (J. Nieves) thanks to DGES contract PB95-1204. 
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Appendix 



The Galilean invariant vertices which appear in the model for eA^ — evrA^, are: 
(a) 'jNN vertex (fig. 3.1(a)): 



^-yNN 



—le 



p + p 



. 3 X g ^Af , 2\ 



(101) 



(b) 7A^A vertex (fig. 3.1(c)): 



Pa 



(S^t X q) 



3 Ma 



(102) 



(c) nNA vertex (fig. 3.1(d)): 

VnNA 

(d) vriVA^ vertex (fig. 3.1(b)): 

fnNN 



I 5'. k ^pa 



(103) 



NN 



B{N, N'n) {ak a {p + p' 



2M, 



(104) 



N 



where q, p, p', Pa J k are the photon, incoming nucleon, outgoing nucleon and pion 
momenta, respectively; ^/s, the invariant energy in the 7* system and M^, and 
Ma are the nucleon, pion and delta resonance masses. In equations ( [103| ) and ( |104D 
we are including the corresponding isospin factors / y B{N, N'tt), respectively. 

Besides the vertices shown in fig. 3.1, in our elementary model for electroproduction 
two more vertices appear: 




(a) (b) 
Fig. A.l 

In fig. A. 1(a) one can see the seagull vertex. It appears from the Ct^nn lagrangian 
via minimal coupling. This vertex is exactly zero for the ir^n and n^p channels and has 
the following expression: 



seagull 



(105) 
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where 



C^'in+n) 



I a{2p + q-k) \ 
2Mn 



\ 



/ a{k - q- 2p) \ 
2Mjv 



V 



C^{7r%) = 



Cin^p) = 



In figure A.l, the vertex (b) corresponds to the ttttj* couphng and it is defined as: 



(106) 



With respect to form factors and couphng constants, their expressions are the 
following: 



We use Sachs Form Factors: 



(107) 



with = 0.71 GeV^] tip = 2.793; /x„ = -1.913. The relationship between Ff (g^) 
(Dirac form factor) and G% is: 



(108) 



4M, 



N 



and = 0. 

For the rest of form factors and coupling constants we take: 



%^ = 0.08 ; ^ = 0.36 ; F^(q2) 



1 



1 - 



(109) 



where Ma = IMGeV. 



-Fy7r7r(? ) 



1-^ 



(110) 



with pI = 0.47^6^2 



- ^ ; A^ ~ 1250 MeV 



(111) 
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where /^(O) = 0.122. It is the 7A^A couphng constant for real photons. 

With respect to the 5* and T operators (transition operator between ^ spin states 

to 2 spin states and respectively between ^ isospin states to ^ isospin states), their 
normalization is: 

<^,M\Si\l,m>^(ll^\XmM) (113) 

< ^,M|ri|^,m >= (l^^AmM) (114) 
where A is a spherical basis index. 
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